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Abstract. We present a general formulation for ray-tracing calculations in curved space-time. The formulation 
takes full account of relativistic effects in the photon transport and the relative motions of the emitters and the 
light-of-sight absorbing material. We apply the formulation to calculate the emission from accretion disks and tori 
around rotating black holes. In our model the emission lines and continuum originate from an accretion disk or 
torus, and the motion of the emitters in the disk/torus is determined by the gravity of the black hole and the space- 
time structure near the black hole. The line-of-sight absorbing medium is comprised of cold absorbing cloudlets. 
These cloudlets are kinematically hot, with their velocity dispersion determined by the local virial temperature. 
The emission from the accretion disk/torus is resonantly absorbed/scattered. Our calculations demonstrate that 
hne-of-sight absorption significantly modifies the profiles of lines from the accretion disks. It is often difficult 
to disentangle absorption effects from other geometrical and kinematics effects, such as the viewing inclination 
and the spin of the black hole. Our calculations also show that emission lines from accretion tori and from thin 
accretion disks differ substantially. Large geometrical obscuration could occur in tori, and as a consequence lines 
from tori generally have much weaker redshift wings at large viewing inclination angles. Moreover, the blue peak 
is truncated. 
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The strong X-rays observed in active galactic nuclei 
(AGN) and some X-ray binaries are believed to be pow- 
ered by accretion of material into black holes. The curved 
space-time around the black hole influences not only the 
accretion hydrodynamics but also the transport of radia- 
tion from the accretion flow. 

Emission lines from thin Keplerian disks around non- 
relativistic stellar objects generally have two symmetric 
peaks (Smak 1969), corresponding to the approaching 
and receding line-of-sight velocities due to disk rotation. 
Because of various relativistic effects, lines from accretion 
disks around black holes do not always have symmetri- 
cal double-peak profiles. The accretion flow near a black 
hole is often close to the speed of light, and emission is 
relativistically boosted. The blue peak of the line there- 
fore becomes stronger and sharper. Moreover, the strong 
gravity near the black-hole event horizon causes time di- 
lation, which shifts the line to lower energies. Emission 
lines from accretion disks around black holes appear to 
be broad, with a very extended red wing and a narrow, 
sharp blue peak (see e.g. the review by Fabian et al. (2000) 
and references therein). Furthermore, gravitational lensing 
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can produce multiple images and self-occultation, further 
modifying the emission line profile. 

Various methods have been used to calculate the pro- 
files of emission lines from accretion disks around black 
holes. The methods can be roughly divided into three cat- 
egories. We now discuss each of them briefly. The first 
method uses a transfer function to map the image of the 
accretion disk onto a sky plane (Cunningham 1975, 1976). 
The accretion disk is assumed to reside in the equatorial 
plane. It is Keplerian and geometrically thin, but optically 
thick. The space-time metric around the black hole is first 
specified, and the energy shift of the emission (photons) 
from each point on the disk surface is then calculated. 
A parametric emissivity law for the disk emission is usu- 
ally used — typically, a simple power-law which decreases 
radially outward. The specific intensity at each point in 
the sky plane is determined from the energy shift and the 
corresponding specific intensity at the disk surface, us- 
ing the Lorentz-invariant property. The transfer-function 
formulation (Cunningham 1975, 1976) has been applied 
to line calculations in settings ranging from thin accre- 
tion rings (e.g. Gerbal & Pelat 1981) and accretion disks 
around Schwarzschild (e.g. Laor 1991) and rotating (Kerr) 
black holes (e.g. Bromley, Chen & Miller 1997). The sec- 
ond method makes use of the impact parameter of photon 
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orbits around Schwarzschild black holes (e.g. Fabian et al. 
1989; Stella 1990; Kojima 1991). The transfer function in 
this method is described in terms of elliptical functions, 
which are derived semi-analytically. The Jacobian of the 
transformation from the accretion disk to sky plane is, 
however, determined numerically via infinitesimal varia- 
tions of the impact parameter (Bao 1992). The method 
can be generalized to the case of rotating black holes by 
using additional constants of motion (Viergutz 1993; Bao, 
Hadrava & Ostgaard 1994; Fanton et al. 1997; Cadez, 
Fanton & Calvani 1998). The third method simply con- 
siders direct integration of the geodesies to determine the 
photon trajectories and energy shifts (Dabrowski et al. 
1997; Pariev & Bromley 1998; Reynolds et al. 1999). 

These calculations have shown how the dynamics of 
the accretion flow around the black hole and the curved 
space-time shape the line profiles. Various other aspects 
of the radiation processes, e.g. reverberation and reflec- 
tion (Reynolds et al. 1999) and disk warping (Cadez et al. 
2003) were also investigated using the methods described 
above. The results obtained from these calculations have 
provided us with a basic framework for interpreting X- 
ray spectroscopic observations, in particular, the peculiar 
broad Fe Ka lines in the spectra of AGN, e.g. MCG-6-30- 
15 (Tanaka et al. 1995). While existing studies have put 
emphasis on the energy shift of the emission, transport 
effects such as extinction have been neglected. Resonant 
absorption (scattering) by ambient material can greatly 
modify the disk emission line profile. This effect was al- 
ready demonstrated in a study by Ruszkowski & Fabian 
(2000), in which a simple rotating disk-corona provides 
the resonant scattering. 

Here, we present ray-tracing calculations of spectra 
from relativistic flows in curved space-time. We include 
line-of-sight extinction and emission explicitly in the for- 
mulation. The radiative-transfer equation is derived from 
the Lorentz-invariant form of the conservation law. It re- 
duces to the standard classical radiative-transfer equation 
in the non-relativistic limit. The formulation can incor- 
porate dynamical and geometric models for the line-of- 
sight absorbing and emitting material. As an illustration, 
we calculate the emision line profiles from thin accretion 
disks and thick accretion tori around rotating black holes. 
The emitted spectra include a power law continuum to- 
gether with a line. This emission is resonantly scattered 
by the line-of-sight-material. We include the contribution 
from higher-order images and allow for self-occultation. 

We organize the paper as follows. In §2 we show the 
derivation of the transfer equation. In §3 we construct the 
equation of motion for free particles in a Kerr space-time 
and for force-constrained particles for some simple para- 
metric models. In §4 we construct a thin disk and a thick 
torus model. In §5 we generalize this by adding in absorp- 
tion due to a distribution of absorbing clouds. In §6 we 
present the results from the models where either emission 
geometry (tori), or absorption (clouds) are important. 



2. Radiative- Transfer Equation 

Throughout this paper, we adopt the usual convention 
c = G = h = 1 for the speed of light, gravitational con- 
stant and Planck constant. The interval in space-time is 
specified by 



dr^ = gafsdx^dx^ 



(1) 



where gap is the metric. 

Consider a bundle of particles which fill a phase-space 
volume element 



dV = dx dy dz dp^ dp^ dp^ 



(2) 



where dx dy dz{= dV) is the three- volume and dp^ dp^ dp^ 
is the momentum range, at a given time t. Liouville's 
Theorem reads 



dV 



(3) 



(see Misner, Thorne & Wheeler 1973), with A here being 
the affine parameter for the central ray in the bundle. The 
volume element dV is thus Lorentz invariant. 

The distribution function for the particles in the bun- 
dle, F(x%p*) is given by 






(4) 



where dN is the number of particles in the three-volume. 
Since dN/dV is Lorentz invariant, F{x'',p^) is Lorentz in- 
variant. From equation ^, we have 



F = 



dN 



p^dVdpdn 



(5) 



where p^ dpdVl — dp'^ dp^ dp^. For massless particles, v — 
c = 1 and \p\ — E. The number of photons in the given 
spatial volume is therefore the number of photons flowing 
through an area dA in a time dt. It follows that 



F 



dN 



E^dAdtdEdn ' 
Recall that the specific intensity of the photons is 
EdN 



" dAdtdEdn ' 

By inspection of equations © and Q , we obtain 



F = 



£3 



(6) 



(7) 



(8) 



where v (= E) is the frequency of the photon. We will use 
this Lorentz invariant intensity, X = i^, in the radiative 
transfer formulation. 

In a linear medium, extinction is proportional to the 
intensity, and the emission is independent of the intensity 
of the incoming radiation. The radiative transfer equation 
is therefore 



as \ V 



(9) 
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where x is the absorption coefRcient, rj is the emission coef- 
ficient and ds is the length element the ray traverses. The 
equation in this form is defined in the observer's frame, 
and the absorption and emission coefficients are related 
to their counterparts in the rest frame with respect to the 
medium via 



e) 



X 



V = [ — 



Xo , 

2 

'70 



(10) 
(11) 



where the subscript "0" denotes quantities in the local rest 
frame. 

The relative energy / frequency shift in a moving 
medium with respect to an observer at infinity is given 

by 

E ^ V ^ '■ ' 



(12) 



-p"Ua\r_ 



(13) 



V ^a |cxD 

where u" is the four-velocity of the medium as measured 
by an observer, and 

ds 
d\ 

The radiative transfer equation (equation [H]) in the co- 
moving frame is therefore 

^ = -p"uM [-Xo(x^ v)I + 770(2;^ v)] (14) 

(see Baschek et al. 1997). 

The results in the co-moving frame can be used to 
determine the intensity and frequency in the other refer- 
ence frames. The ray is specified by choosing x"(Ao) and 
p"(Ao). From the geodesic equation, we have dp^/dX + 
^'a^P P^ = 0, where we have scaled A by m for massive, 
and by 1 for massless particles. The total derivative of I 
is therefore 



dA 



dl dx° 



dl dp° 



dx" dX dp" dX 



dl 



dl 



P dx- ^^''PP dp- 



(15) 



This, combined with equation H14(l . yields 



= P 



, dX 
dx" 



i/37-p p 



dl 
dp^ 



(16) 



which is the same as that derived by Lindquist (1966) from 
the Boltzmann Equation. 

The metric and the initial conditions define the rays 
(the photon trajectories in 3D space), and the solution can 
be obtained by direct integration along the ray. For sim- 
plicity, we assume the refractive index n = 1 throughout 
the medium. The solution to equation (|14() is then 

X(A) = X(Ao) cxp ( /" xo(A',i^o)«aP"dA'J (17) 

cxp / xoi>^",'^o)uaP°'dX" ]T]o{X',iyo)uaP°'dX' . 



In the non-relativistic limit, Uajf^ = 1, and the equation 
recovers the conventional form (see Rybicki & Lightman 
1979). 

3. Particle Trajectories 

3.1. free particles 

To determines the photon trajectories we need to spec- 
ify the metric of the space-time. We consider the Boyer- 
Lindquist coordinates: 



dr' = l- 



2Mr 



dr 



AaMrsm^e , ,, S , o 
dtdd —dr 

E ^ A 



-^d0^ - 



2o?Mr sin^ 9 



sin^ edct)"^ , (18) 



where M is the black hole mass, with the three vector 
(r, 9, (f>) corresponding to spherical polar coordinates, and 
defining T, = r'^+a'^ cos^ 9 and A = r^ - 2Mr + a'^. The di- 
mensionless parameter a/M specifies the spin of the black 
hole, with a/M = corresponding to a Schwarzschild 
(non-rotating) black hole and a/M = 1 to a maximally 
rotating Kerr black hole. 

The motion of a free particle is described by the 
Lagrangian: 



C 



2Mr 



AaMr sin"' ■ 



-i(j) 



-Er 



2a^ Mr sin^ 



sm 



A^ 
2 /) 12 



(19) 



(here x" — dx"/dX). The Lagrangian does not depend 
explicitly on the t and (f) coordinates. The momenta in the 
four coordinates are therefore 



Pt - dC/di - 

Pr = dC/dr = 

pe = dC/d9 = 
p^ = dC/dcj) = 



-E , 



A 
Y.9 
L . 



(20) 

(21) 

(22) 
(23) 



with E being the energy of the particle at infinity and 
L the angular momentum in the direction. The corre- 
sponding equations of motion are 



t 



E 



2r(r2 



)E - 2aL 



hi 



2arE- 



EA 
Ei- L^- S^2 

(E-2r)i/sin 



H)a^ cos^ 9- L^ cot^ i 



EA 



(24) 
(25) 
(26) 
(27) 



where Q is Carter's constant (Carter 1968), and H is the 
Hamiltonian, which equals for photons and massless par- 
ticles and equals —1 for particles with a non-zero mass. 
(See Reynolds et al. (1999) for more details.) For simplic- 
ity, we have set the black-hole rest mass equal to unity 
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(M = 1) in the equations above. This is equivalent to 
normahzing the length to the gravitational radius of the 
black hole (i.e., set Rg = GM/c^ — 1), and we will adopt 
this normalization hereafter. 

There are square terms of f and 6 in two equations 
of motion. They could cause problems when determining 
the turning points at which f and 6 change sign in the 
numerical calculations. To overcome this, we consider the 
second derivatives of r and 9 instead. From the Euler- 
Lagrange equation, we obtain 



"= E 



Ari]-2r2.2 (r-l)E-rA 



E2 



-t^ + 



sin^ e 



-2a sin 



= — < sm u cos ( 



A2 
S-2r2 

E2 

-2r2 - 



rr 



a^ SIT? ( 



S2 



2a^r ■ 



-t' - 



-t(j)- 
Aar{r 



2a? sin 6 cos 



+a^ 



A + 2r(r2 + a^f 



-t^^-r^ 



2rfe 



(29) 



In terms of the momenta and the Hamiltonian, the equa- 
tions above can be expressed as 



+2r{r'^ + a'^)E^ - 2aEL] 



l-rAiJ 
2p,2(r-l) 



sm f cos 



sm 



-a\E^ + H) 



(30) 
(31) 



where k = Q + L^ + a'^{E^ + H). Equations (|1TJ, if^ . 
(|23I) , l|23 , EOl and (EU are the equations of motion. 

3.2. Motion in the presence of external forces 

The equations of motion obtained in the previous section 
are applicable to free particles only. In a general situation 
external (non-gravitational) forces may be present and we 
need to specify the external force explicitly in deriving the 
equations of motion. However, in the setting of accretion 
disks around black holes we can often treat the effect of 
the external force implicitly which we will discuss in more 
detail in the following subsections. 

3.2.1. Rotational and Pressure Supported Model 
Here we consider a simple model such that 

i></)>r>^. (32) 

As r and 9 are small in comparison with other quantities, 
they can be neglected as a first approximation. 
The equation of motion reads 






(33) 



where a^ is the four- acceleration due to an external force 
per unit mass. For axisymmetry (which is a sensible as- 
sumption for accretion onto rotating black holes), d/dt/i = 
and a? ~ 0. The identity u°'aa — together with 
r = and (? = imply that a* = 0. We may also set 
a*" — for simplicity. Because we have an extra equation 
from the identity m"uq = —1, a^ can be determined self- 
consistently under the approximation = 0. This scenario 
thus corresponds to flows supported by rotation in the f 
direction and by pressure in the 9 direction. 

Inserting the affine connection coefficients for the Kerr 
metric into the equation of motion yields quantities iden- 
tical to zero on the left hand side of the equations for the t 
and directions. This leaves only the non-trivial momen- 
tum equation in the radial direction: 



(28) = - 



I]-2r2 

Y? 



r + 2 



- I r + a^ sin^ ( 



I]-2r2 

E2 

S-2r2 



S2 



sm 



which further simplifies to 



rE2 sin^ i 



^^ - [i- a sin^ I 



2r2 - S ' 
Here we choose the positive solution 



/ \/rS sin 9 . 9 „ 



V\/2r2-E 



(34) 



(35) 



(36) 



which corresponds to the same rotation as the black hole. 
This solution thus allows the flow to match the rotation 
of a prograde accretion disk. 
From the metric we have 

2r \ ., 4ar sin^ 9 ■ ■ 



r"^ +a^ + 



2a^rsii?9 



,2 q12 



(37) 



Combining equations (|36|l and H37|) yields 



2r2 



^/2r 



(38) 



It follows that the components of the four-velocity are 

i = i (y.^ + as\xv9^/2r'^ - y\ , 

f = , 
^ = 0, 



v/272" 



Csin6' 



(39) 



where 



C = WS (E(r + 1)- 4r2 + 2asm9^Jr{2r'^ -Y)\ . (40) 

The marginally stable orbit for particles is defined by 
the surface where 

dE , , 

= . (41) 



dr 
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From equations fl^i . ifTTI) and (|^ . we have 
E^-({j:~2r)^ + asmey/2r^-j:\ . (42) 



c 

After differentiation, we remove the non-zero factors in 
the expression and obtain the condition 



AS^ - 4r (2r2 - S) ( ^/2r^ - S - a sin 9y/^) =0. (43) 



Setting a = gives r — 6, which is often regarded as the 
hmit for the inner boundary of an accretion disk around 
a Schwarzschild black hole. This value is the same as 
that derived by Bardeen, Press & Teukolsky (1972) us- 
ing d'^p'^/dr^ = 0. 

Before we proceed further, we must note that the ex- 
pressions for the velocity components in equations (|39|l 
hold only for regions "sufhciently" far from the black-hole 
event horizon. The approximation that we adopt in the 
model breaks down when the square root in the denomi- 
nator approaches zero. This occurs at the light circulari- 
sation radius rdr, which is given by C = 0, or equivalently 



S(r + 1) - 4r2 + 2asm9y/r{2r'^ - S) 



= 



(44) 



Moreover, the assumption of ^ = is also invalid for radii 
smaller than the radius of the marginally stable orbit — 
the flow is neither rotational nor pressure supported and 
it follows a geodesic into the event horizon. 

3.2.2. Isobaric Surfaces 

In a stationary accretion flow, the acceleration must be 
balanced by some forces, e.g. the gradient of gas or ra- 
diation pressure. As the local acceleration a" can be cal- 
culated from the rotation law uj{r, 9) we can derive a set 
of isobaric surfaces when a rotation law is given. For a 
barotropic equation of state of the accreting matter the 
isobaric surfaces coincide with the isopicnic (constant- 
density) surfaces. 

The accelerations in the f and 9 directions are 



E , ^-2r\ ^ . 

—a = 7^ — |r — asmf 

A S^ 



r sm 



2 /) 12 



-Ea = sin 6' cos ( 



2r 

£2 



ai - (r^ + a^ 



A02 



The surface of constant acceleration is given by 



dX 



(45) 
(46) 



(47) 



(here and hereafter dxf^^^.f/dX = dx" /dX\x^^^i)- The sta- 
tionary condition implies dt/dX — 0, and axisymmetry im- 
plies d(t)/dX — 0. Without losing generality, we can choose 
t = (f) = on the surface. Thus, equation (|47|l becomes 



= 



Ea'' drsuTf 



A 



dX 



dX 



-P2 



d9. 



,d9 



surf 



dX 



dX 



(48) 



where 



(3i = 



E-2H 



E2 
P2 = sin cos 



a sm ( 



r sm 



2r /a , 2 , 2 ' ^ 



A+^ --(^^+«^) 



(49) 



and drsuTf/dX and d9suif/dX determine the intersection 
of the isobaric surfaces and the (r, 9) plane. By rescal- 
ing equation (|48|) with a factor of -^/A/E and making use 
of the invariance 



'df 
we obtain 

drsurf 



Er 



(50) 



01 



dX' ^pl + A/32 



d9^ 



-P2 



dX' 



VMTW, 



(51) 



These two differential equations can be solved numer- 
ically and yield the isobaric surface as a parametric func- 
tion of A'. 

4. Model Accretion Disks and Tori 

We now demonstrate using the equations of motion above 
to construct the emitter models. The first is a geomet- 
rically thin accretion disk, in which the emitting parti- 
cles are in Keplerian motion. The second is a torus, a 
3-dimensional object with non-negligible thickness. 

4.1. Accretion Disk 

When space-time curvature is important, the Keplerian 
angular velocity of a test particle around a gravitating 
object is no longer ojk — r~^"^ ^i^q expression in flat 
space-time. Instead, the Keplerian angular velocity in a 
plane containing the gravitating object can be obtained 
by setting 6 = tt/2 in equations (|39|l and H4()|l . Hence, the 
components of the four- velocity of the particles in the disk 
are 

t = 



T\l r"^ — 3r + 2a^ 



f = , 
^ = , 



1 



3r + 2a^ 



(52) 



and the rotational velocity of a Keplerian accretion disk 
around a black hole is 
1 



Wfc 



r3/2 



(53) 



(Bardeen, Press & Teukolsky 1972). 

The relative energy shift of the emission between the 
disk particle and an observer at a large distance is de- 
termined by equation (|12|l . with u" as given in equation 
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(|52l) . Keplerian disk images can be found in many exist- 
ing works (e.g. Bromley, Miller & Pariev 1998), and we 
do not show disk images here. The general characteristics 
are that a disk image is asymmetric, with the separatrix 
for the energy shift of the emission no longer bisecting 
the disk image into two equal sectors, one for red shift 
and another for blue shift. The whole disk appears to be 
reddened, especially at the inner rim. 

4.2. Accretion Torus 

To determine the geometry and structure of an accretion 
torus self-consistently is beyond the scope of this paper. 
Here, we consider a simple parametric model, with an an- 
gular velocity profile given by 



1 



(rsin6')3/2 + a 



V r sin 9 ) 



(54) 



The quantity r\ is the radius (on the equatorial plane) at 
which the material moves with a Keplerian velocity. The 
parameter n adjusts the force term, such as a pressure 
gradient, to keep the disk particles in their orbits, and it 
determines the thickness of the torus. In this study we 
just take n = 0.21 without losing generality. If the torus 
is supported by radiation pressure, its inner edge is de- 
termined by the intersection of the isobaric surface with 
either one of two surfaces. These two surfaces provide the 
constraints, inside which the pressure-supported solution 
does not hold. The first is a surface defined by the orbits 
of marginal stability. For w(r, 0), it is given by 



== 2asin^ 



sm 



6r(r^ 



2r 
1 

S 

-|-Asm^0 u- 

or 



u?- 



a^r sin^ ( 



-f 3A-S 
6ar sin^ Q 



E-2r^ 
-2r2 



S-2r2 



1-^ 



S-2H 

Y? 



(55) 



The second is the limiting surface where the linear velocity 
approaches the speed of light. It is given by 







S - 2r -f 4arw sin^ Q 



((^ 



)I] + 2a2rsin2 



Lo sin" ( 



(56) 



Usually the former is larger than the latter. The outermost 
of these two surfaces determines the inner boundary and 
hence the critical surface of the torus. 

Figure^shows the critical density surfaces of two tori. 
The first torus is around a Schwarzschild black hole and 
the second torus is around a maximally rotating black 
hole. The tori are constructed such that their specific an- 
gular momentum has a profile similar to those of the simu- 
lated accretion disks in Fig. 3. of Hawley & Balbus (2002). 

In our model the boundary surface of the torus is de- 
termined by a single parameter, n, which specifies the in- 
dex of the angular-velocity power law. Its value is selected 



such that the angular- velocity profile matches the profiles 
obtained by the numerical simulations — here we consider 
that of Hawley & Balbus (2002). Model tori can be con- 
structed using various different methods. An example is 
that in a study of dynamical stability of tori around a 
Schwarzschild black hole carried out by Kojima (1986), 
the model parametrizes the angular momentum instead 
of the angular velocity. We note that the aspect ratios of 
the torus surfaces obtained by Kojima (1986) and those 
shown in Fig. ^ are similar. 
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Fig. 1. The critical emitting surface of the torus (solid 
line), the surface defined by the orbits of marginal stabil- 
ity (dashed line) and the limiting surface where the linear 
velocity derived from equation H54|) approaches the speed 
of light (dotted line) for a Schwarzschild black hole (top 
panel) and a Kerr black hole with a = 0.998 (bottom 
panel). In both cases we consider an profile index n = 0.21 
(see Equation |54jl. which gives the tori an angular mo- 
mentum profile similar to those obtained in the accretion 
disk simulation of Hawley & Balbus (2002). The Kepler 
radius r-k is 15 for the Schwarzschild black hole and is 12 
for the Kerr black hole. 



5. Extinction 

The generic setting of the system under our investigation 
is that emitters with various strengths are distributed in 
space in a curved space-time, and the radiation is attenu- 
ated, and may be re-emitted, when propagating. The emit- 
ters and the line-of-sight material are in relativistic motion 
with respect to the observer and also with respect to each 
other. An example is that shown in Fig. [5], in which the 
emitters are the surface elements of an accretion disk and 
the absorbers are some clouds in the vicinity of the disk. 
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The photon trajectories and the motion of the emitters 
and absorbers are affected by the space-time distorted by 
the central black hole. 

To construct the model we need to determine 

(i) the rays that connect the emitters, absorbers and ob- 
server, 

(ii) the four- velocities of the emitters and absorbers /scat- 
tcrers, 

(iii) the spatial distributions of the emitters and the ab- 
sorbers/scatterers, and 

(iv) the effective cross section of the absorbers/scatterers. 

In the previous section, we have shown how to obtain (i) 
and (ii); in this section we incorporate (iii) and (iv) into 
the radiative-transfer calculations. 



Distribution of Absorbing Clouds 

,-'„« O '"■-, 
/' O Qf^lo \ 




Black Hole 
Accretion Disk 

Fig. 2. A schematic illustration of the physical setting 
on which the radiative-transfer formulation is constructed. 
The space-time curvature is determined by a gravitating 
body, which is a rotating black hole. The emitting ob- 
ject is the accretion disk around the hole. The cloudlets 
in the vicinity attenuate the radiation from the accretion 
disk. The properties of the radiation are specified by the 
4- velocity of the photon, m" The other three important 
quantities in the model are the 4-velocity of emitting ele- 
ments in the accretion disk u^j^j, the bulk velocity of the 
absorbing cloudlets u"[ and the "microscopic" 4-velocities 
of the cloudlets u'^^ with respect to the bulk motion of the 
cloudlets. 



5.1. An illustrative model 

We consider a model with the geometry shown in Fig. [21 
The photons are emitted from the elements on the top and 
bottom surfaces of a geometrically thin disk in a Keplerian 
rotation around a Kerr black hole. The radiation is reso- 
nantly scattered (absorbed) by plasma clouds and is at- 
tenuated in its propagation. We omit emission from the 
clouds (i.e. ryo = in equation ^3]) in this study. The 
size of the clouds is small in comparison with the length 
scale of the system. They are not confined to be in the 
equatorial plane and are in orbital motion, supported by 
some implicit forces (which may be radiation, kinematic 
or magnetic pressure gradients). These clouds have a large 



(thermal) distribution of velocities, in addition to their 
collective bulk velocity. 

We assume that the radiation scattered into the en- 
ergies of the lines is insignificant and ignore the photons 
that are scattered into the line-of-sight. Under this approx- 
imation, scattering simply removes the line photons and 
causes extinction similar to true absorption. Thus, for sim- 
plicity, hereafter we do not distinguish between scattering 
and absorption, ^ and the two terms are interchangeable, 
unless otherwise stated explicitly. 

The rays originating from the accretion-disk surface 
that can reach the observer are determined by the 4- 
momenta of the photons, which are calculated using equa- 
tions (Hill J (113 and (|2nil- The 4-velocities of the emitting 
surface elements on the accretion disk are given by equa- 
tions H52|) . These determine the relative energy shifts of 
the photons between the emitters and the observer. What 
we need next is to determine the relative energy shifts 
between the emitters and the absorbers. Then, we need 
a model mechanism by which the absorption takes place, 
and to derive the resonant absorption condition for the 
absorption coefficient. 

Now we construct a model for the spatial distribution 
and the velocities of the absorbers. Consider a paramet- 
ric model in which the bulk 4-velocities of the clouds are 
given by equations (|39|l and l|40|l . In this model the bulk 
velocities of the clouds in the equatorial plane matches the 
4-velocities of the accretion disk. 

The clouds themselves are cold, and the thermal ve- 
locity of the gas particles inside are much smaller than 
their bulk motion and root-mean-square velocity disper- 
sion. However, the clouds have a large velocity dispersion, 
given by the local virial temperature, which is comparable 
to the energy of the emission lines of interest. Therefore, 
the clouds can be considered as relativistic particles in the 
calculation. Using the bulk-motion velocities obtained by 
(|39|l and H40I) together with the virial theorem, we can 
derive this temperature and determine the velocity distri- 
bution of the absorbing clouds. 

The clouds fill most of space, with a radially depen- 
dent number density. However, close to the black hole, 
the assumptions above break down, and the axial force 
cannot support the clouds out of the equatorial plane. 
When this happens, they will flow along geodesies directly 
into the hole. In the numerical calculation, we determine 
the asymptotic boundaries at which the left hand sides of 
equations l|ll|l and (|1^ vanish. This is done by evaluating 
these expressions and testing to see if they are negative 
at each point along the photon rays. Inside that surface, 
the number density of the clouds will be much less than 
that outside, which we approximate by setting it to zero 
in this zone. 



Here, in spite of the term "resonant scattering" (which de- 
scribes the nature of the process), we have ignored the scat- 
tering of continuum photons into the line energy bins and the 
energy (frequency) redistribution in our calculations. 
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5.2. The absorption coefficient 

We assume that the absorption is due to "cold" cloudlets 
with high virial velocities. The absorption coefficient of 
individual cloudlets is 



Xi oc cr 



U^kg +£^linc 



(57) 



(with a as the effective absorption cross section of the 
cloudlet, and ka the photon four-momentum). The ab- 
sorption rest frequency is E'line, and E^ is the energy of 
the photon in the bulk rest frame. The total effective ab- 
sorption coefficient xo is the sum of the contribution of 
these cloudlets, i.e.. 



XO 



E 



Xi 



Converting the sum into an integral in momentum space 
yields the absorption per unit length in the rest frame as 



-27rAo- 



p dpd^exp{~~E/Q)u°'kaS 



E~t 



, (59) 



where we have defined /i = cos 0, A is a normalization 
constant, Q is the temperature in relativistic units (with 
ks = 1), and E and p are the energy and momentum of 
a gas particle in the bulk rest frame. This form assumes 
isotropic thermal motion in the rest frame. 

There are three terms that we need to determine be- 
fore we can evaluate the integral, and they are the nor- 
malisation parameter A, the temperature Q and the pho- 
ton energy in the rest frame of the absorbing particle 
u°'ka. When these variables are determined, we can then 
parametrise a after integration is carried out. 

5.2.1. The normalisation parameter A 
We now derive A starting from 



/>oo 

N = Att\ l p^dpexp i-E/Q) 
Jo 



(60) 



where A'^ is the number density of absorbing clouds. 
Integrating yields 



A = 



TV— 



(61) 



where K^{x) is a modified Bessel function, and m is the 
average cloud mass. This is called the Jiittner distribution 
and corresponds to Maxwell's distribution except in the 
case of a relativistically high temperature. 



/•OO 

Etot = 4^A / p'dpE exp {-E/Q) 
Jo 



(62) 



Integrating this yields the energy per unit mass as 



Eu 



^3(t) e 



(63) 



Nm i^2(f) 

Using conservation of energy and angular momentum, 
we calculate the thermal energy of the virialised relativis- 
tic gas of absorbing clouds. 

At infinity the medium has 

Einit = Nm , (64) 

iinit = ifin • (65) 

Close to the black hole it has 

Nm r 

-C/firi 



(58) Lfin = 
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(E - 2r) Vr + a sin 9y/2r'^ - E 
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(66) 
(67) 



The energy released by the gas falling from infinity 
and slamming into a wall moving with a velocity given by 
equation Ip!?^ is 

— Etot — u'^Pa — — -Efin^init " ifin'/'init • (68) 

After simplification, this becomes 
^tot 1 



Nm C^ 



{2r^ - E)(r2 + a^) - 2ar sin 9 y/? \/ 2r^ ~ T, 



(a sin6l\/2r2-E + T,y/A ( 



(69) 



Thus the temperature of the media can be derived 
using equations H63(l and H69() . Unfortunately, this yields 
an implicit relation of to/0 that contains transcendental 
functions. The modified Bessel functions can be expanded 
in the limit where O <C to, which corresponds to an "al- 
most relativistic" gas. Since the potential energy released 
in accretion is of the order of a few percent of the rest 
mass of the infalling material, this approximation should 
hold in AGN. 

Expanding to second order in 8/to, cancelling the ex- 
ponential factors, and then solving the resulting quadratic 
yields 




10 fEt, 



3 \Nm 
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(70) 



Thus we have an explicit description of how the kinematic 
temperature varies with position. 

5.2.3. The photon energy u'^ka 

In the rest frame, the motion of the thermalised medium 
is isotropic. Thus we can simplify the problem by aligning 
an axis along the photon propagation vector and working 
in a local Lorentz frame so that 



5.2.2. The temperature G 

The total energy in the distribution of clouds is given by Thus 



A:„==£;^(-l, 1,0,0) , 

and 

P" ^mu" = {E,pii,py,p^) 
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(71) 
(72) 

(73) 
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5.2.4. Evaluation of the 5— function 



Using the relation that 
d 



P_ 
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(74) 



we obtain 
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We change the variable pdp to EdE. The integral in equa- 
tion (|59|l can now be simplified to 
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(76) 



Integrating this gives the absorption coefficient in the rest 
frame of the gas 
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This equation can be recast in terms of x by using equa- 
tions (dni) and IpH) . 

5.2.5. The effective absorption Na 

The absorption coefRcient depends upon N, the number 
density of the clouds, and on a, the absorption cross- 
section per cloud. These are in general a function of po- 
sition. Since the pressure gradient and the inflow velocity 
in the f and directions are ignored in our approximation 
of the flow, the density profile cannot be determined in a 
fully self-consistent manner via the mass continuity equa- 
tion. To overcome this, we assume a simple two-parameter 
profile 



Na — fj^r 



-P 



Where ctq can be considered as a proportionality constant 
fixing the density and opacity scales. Without losing gen- 
erality, we adopt a value (3 = 3/2. The results can then 
easily be generalized to other values of (3. For the case 
of a swarm of absorbing clouds around a black hole, the 
optical depth is approximately one. This corresponds to 
the case where ctq is of order 0.5 or greater, with the in- 
tegration proceeding radially to the event horizon. The 
effective optical depth depends greatly on the paths the 
photons take. 




Fig. 3. The profiles of emission lines from thin Keplerian 
accretion disks around a Schwarzschild black hole (top) 
and a Kerr black hole with a ~ 0.998 (bottom). The line 
profiles are normalized such that the flux F{E) = 1 at 
E/Eq = 1. The viewing inclination angle is 45°. The in- 
ner radius of the accretion disk is at the last stable orbit, 
and the outer radius is 10i?g. The line emissivity on the 
disk surface is a power-law which decreases radially out- 
ward, and the powerlaw index is —3. Only emission con- 
tributed by the direct disk image is considered; emission 
from higher-order images is not included. The emission 
is unabsorbed. In each case, the solid lines correspond to 
the line profiles obtained by semi-analytic calculations de- 
scribed in Fanton et al. (1997), and the circles represent 
the line profile obtained by our numerical ray-tracing cal- 
culations. 



(78) 6. Spectral Calculations 



6.1. No absorption 

6.1.1. Accretion disks 

We calculate the observed energy of the flux from a point 
on the planar disk using equation H12|l . We ignore absorp- 
tion. The line emissivity has a power law profile, which 
decreases radially from disk centre. The intensity is pro- 
portional to the third power of the relative shift; the flux 
has an extra factor of vq/v due to time dilation, i.e. it is 
proportional to the fourth power of the relative frequency 
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Fig. 4. Comparison of the line profiles for cases considering only the direct disk image (dotted line) and cases including 
high-order images (solid line). The line profiles are normalized such that F{E) = 1 at E/Eq = 1. The viewing inclination 
angles of the disks are 75°, 85° and 89° (panels from top to bottom). Line profiles for Schwarzschild black holes are 
shown in the left column; line profiles for the Kerr black holes with a = 0.998 are shown in the right column. The 
outer radius of the accretion disk is 20i?g, and the index of the cmissivity powerlaw is —2. 



shift. (Note that F in §2 is the distribution function, not 
the flux of the emission.) 

Our calculation reproduces the line profiles of direct 
images of accretion disks as those shown in Fabian et 
al (1989), Kojima (1991), Fanton et al (1997), Bromley 
et al (1997), and Reynolds et al (1999). Figure shows 
two example line spectra calculated using the method de- 
scribed above. The spectra contain only emission from the 
direct image. We also show line spectra obtained using the 
method by Fanton et al (1997) for comparison. The results 
in the two ca Iculations arc in excellent argeement. 



We also carry out spectral line calculations which in- 
clude contribution from higher-order disk images. (Here 
and hereafter we assume that the emissivity powerlaw has 
an index of —2, except where otherwise stated explicitly, 
i.e. I(Ao) oc r^^) Our calculations show that the con- 
tribution of the higher-order images are significant only 
at high inclination angles (see Fig.Q. The emission from 
high-order images is mostly at frequencies close to the rest 
frequency of the line, because the region where highly red 
and blue-shifted emission originates is obscured. 
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Fig. 5. Images of an accretion torus around a Kerr black 
hole viewed at inclination angles of 45° (top) and 85° (bot- 
tom). The angular momentum of the black hole a = 0.998. 
The torus is constructed using the velocity law in §4.2. The 
inner and outer radii are determined by the critical sur- 
face as described in the same section. High-order images 
are shown. The energy shift of the emission is represented 
by a false-colour map laid on the torus surface. Red-shifted 
emission is coloured yellow/red and blue-shifted emission 
is coloured blue. The grey lines indicated by the arrows 
mark the region where the emission has zero energy shift, 
separating the torus into regions with red and blue shifted 
emission. 

6.1.2. Rotational Torus 

We now investigate the emission from an accretion torus. 
We consider a model in which the inner radius of the torus 
is determined by the marginally stable orbits of the par- 
ticles. These marginally stable orbits, which depend on oj 
and duj / dr, form a surface in a three-dimensional space. 
The marginally stable orbit for particles in Keplerian mo- 
tion in the equatorial plane is 6 R^ around a Schwarzschild 
black hole and is 1.23 i?g around a Kerr black hole with 
a = 0.998. In equation (|54|l . Keplerian motion corresponds 






0.5 



1 
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Fig. 6. Profiles of emission lines from an accretion torus 
around a Kerr black hole with a = 0.998 at viewing in- 
clination angles of 45° (solid line) and 85° (dotted line). 
The parameters of the torus are n — 0.21 and r^ = 12i?g 
in the velocity profile equation (|54|l . 




Fig. 7. A comparison of the profiles of lines from an accre- 
tion disk (dotted line) and an accretion torus (solid line) 
around a Kerr black hole with a = 0.998. The disk and 
the torus are viewed at an inclination angle of 85°. The 
parameters of the torus are the same as those in Fig. El 
The disk has an outer radius of 20i?g, an inner radius of 
1.23i?g, and a Keplerian velocity profile. 



to the case with the index n = 0. In the torus model 
that we consider here, n is not necessarily zero. As a con- 
sequence, the location of the marginally stable orbit of 
particles in the tori and conventional Keplerian disks are 
different. (See Fig. Q] in §4.2). We use a surface-finding 
algorithm to determine the boundary of the torus (see 
Appendix 0. 

In Fig. El we show three-dimensional images of the 
model torus around a Kerr black hole. We include the 
first four image orders. The inclusion of high order im- 
ages is mandatory due to the 'mixing' caused by the ex- 
tension out of the equatorial plane (Viergutz 1993). The 
torus is viewing inclination angles of 45° and 85° (top and 
bottom panels respectively). The left-right asymmetry is 
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caused by inertial-frame dragging. The multiple images 
are consequences of gravitational lensing. At small incli- 
nation angles, only the surface above the equatorial plane 
of the torus is seen in the direct image. At very large in- 
clination angles, the surface below the equatorial plane is 
severely lensed and also becomes visible. 

The false-colour map laid on the torus surface show 
the energy shift of the emitted photons (determined by 
equation 22 )j ^s viewed by a distant observer. The sepa- 
ratrix, which corresponds to zero energy shift, divides the 
torus surface into regions of blue energy shift and regions 
of red energy shift. At large inclination angles the inner 
surface of the near side of the torus is not visible, and the 
inner surface of the far side is obscured by the near side 
of the torus. Thus, emission with the largest energy shifts 
is hidden. This is very different to the situation for a pla- 
nar accretion disk - regardless of the viewing inclination 
and the visual distortion, the emission from the innermost 
part of the disk is always visible. 

Figure El shows the resulting line profiles obtained by 
integrating the emission over the images shown in Fig. [S] 
The line profile of a torus viewed at 45° is similar to that 
of the planar accretion disk. It has a sharp blue peak and 
smaller red peak. It also has an extended red wing. This is 
due to the fact that the projections of a torus and a disk 
on the sky plane are very similar at low inclination angles. 

However, our calculations show that geometric effects 
are very important for large viewing inclination angles. 
When some part of the emission region is self-obscured, 
the resulting profile, as observed from infinity, to be com- 
pletely different from that of the fiat disk (see Fig.|3). For 
a torus with large viewing inclination angles, the inner sur- 
face of the torus tends to be obscured. This corresponds 
to the region where the most redshifted flux of the line 
is emitted (due to large transverse red shift and gravita- 
tional red shift). This makes the red wing less prominent. 
The outer surface of the torus is visible from all inclina- 
tions. As a result the line profile tends to be singly peaked, 
with the maximum at the unshifted line frequency due to 
the emission from the outer surface dominating. By alter- 
ing the geometry of the emitter a wide variety of emission 
profiles can be obtained. 

6.2. Resonant Scattering 

We use a disk model to illustrate the resonant scattering 
effects. We assume the inner edge of the accretion disk 
is given by the marginally stable orbit. We use an outer 
disk radius of 20i?g in all the disk simulations. This was 
chosen to accentuate the relativistic effects. The emission 
line profile is assumed to be a delta function. We collate 
the light from the first four image orders of the accretion 
disk. 

The line profiles are obtained from 750 x 750 pixel 
images. The intensity scale on the graphs is in arbitrary 
units. (It is just the log of the sum of pixel intensity over 
the image, as a function of frequency.) The x-axis of the 



graph is in units of the line rest energy, with E/Eq = 1 
corresponding to the unshifted line. 

We bin the injection spectrum and the absorption co- 
efficient linearly with energy. There are 1000 bins from 
E = to E — 2Eq. Since the thermal width of the line 
is small, we assume that the emission line is narrow, with 
a width of only one energy bin. After the calculation, we 
smooth over sets of ten adjacent bins to remove numerical 
noise which has a large scatter. See Appendix for more 
information about the algorithm used. 

We investigated two space-time models. One with a ~ 
corresponding to a Schwarzschild black hole, and one 
with a — 0.998, corresponding to a maximally spinning 
Kerr black hole. We have plotted spectra containing the 
continuum and the continuum plus line as we have varied 
the opacity of the absorbing clouds. 

We have also shown how the spectra change with incli- 
nation due to geometric effects. We have included a nearly 
edge-on model (i = 85°), and a model with a moderate 
inclination of 45°. See Fig. |H1 for the results of absorbing 
the power law continuum, and Fig. El for the results of 
absorbing both the line and continuum. 

7. Discussion 

We have modelled line profiles from accretion disks to 
demonstrate the use of a general formulation for transfer 
of radiation through relativistic media in arbitrary space- 
times. In this paper, we used the transfer of emission from 
AGN as an illustration. In this model, we parametrized the 
disk/torus to describe the emitters and the space-density 
distribution of the absorbers. We took into account rela- 
tivistic effects on the bulk dynamics and the microscopic 
kinematic properties of the absorbing medium. 

Resonant absorption/scattering of line emission from 
accreting black holes in the general relativistic framework 
had been investigated previously by Ruszkowski & Fabian 
(2000). In their study, a thin Keplerian disk was assumed 
and the absorbing medium is a spherical corona of con- 
stant density centred on the black hole. The corona is ro- 
tating, with local rates obtained by linear interpolation 
from the rotation rate of a planar Keplerian accretion 
disk and the rotation rate at the polar region caused by 
frame dragging due to the Kerr black hole. The Sobolov 
approximation was used in the resonant absorption cal- 
culations, and a Monte Carlo method determined the re- 
emission/scattering. 

Our calculation is different to that of Ruszkowski & 
Fabian (2000) in the following ways. Firstly, the emitters 
are not confined to the equatorial plane, i.e. they can be 
thin accretion disks or thick tori. Secondly, the absorbing 
medium is a collection of (cold) clouds with relativistic 
motions. The number density distribution of the clouds is 
parametrized by a powerlaw decreasing radially. The local 
bulk (rotational) velocity of the clouds is determined by 
general relativistic dynamics, and the velocity dispersion 
is calculated from the Virial theorem. Thirdly, we do not 
assume the Sobolov approximation. The resonant condi- 
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Fig. 8. Absorbed continuum emission from accretion disks around Schwarzschild and Kerr black holes with a = 0.998 
(left and right columns respectively), for i = 45° and 85° (top and bottom panels respectively), uo = 0, 0.05, 0.2 and 
0.5 (curves from top to bottom in each panel). 



tion for the absorption coefficient is derived directly from 
the kinematics of the absorbing cloud particles. Fourthly, 
we ignore the contribution from re-emission to the line 
flux. However, we include emission from higher order disk 
images, in addition to the direct image. 

One of the main differences between the two studies 
is the treatment of resonant absorption. In the Sobolov 
approximation, the absorption takes place locally (see 
Rybicki and Lightman 1979). The line profile is practically 
a delta function; otherwise, the assumption of quasi-local 
absorption breaks down. Moreover, it requires that the 
absorbing medium is a radial flow. The emission lines are 
broadened because of relativistic effects, and the motion 
of the medium is rotationally dominated. The locality of 
the absorption (required by the Sobolov approximation) 
therefore breaks down. To overcome these difficulties, we 
abandon the Sobolov approximation but, instead, employ 
full ray tracing. 



In modelling the bulk flow, we derived the equation of 
motion of the absorbers using the rotational-support ap- 
proximation for the accretion disk. In addition, we assume 
that the flow is supported out of the equatorial plane. 
There is negligible radial force in the equations of motion, 
and the radial velocity can be neglected. To include out- 
flowing as a wind, or inflowing, requires additional com- 
ponents in the equation of motion. This complicates the 
formulation, in particular, when matching the flow bound- 
ary condition at the surfaces of the accretion disk/torus. 
A self-consistent boundary condition requires a dissipa- 
tion mechanism in the boundary regions of the disk/torus 
and absorbing clouds. The inclusion of such dissipation 
is beyond the scope of this paper and this issue will be 
addressed in future works. 

As the rays propagate from the accretion disk to the 
observer, they experience position-dependent absorption. 
The absorption depends upon the velocity profile of the 
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Fig. 9. Line and continuum emission from accretion disks around Schwarzschild and Kerr black holes with a = 0.998 
(left and right columns respectively), for i = 45° and 85° (top and bottom panels respectively), ctq = 0, 0.05, 0.2 and 
0.5 (curves from top to bottom in each panel). 



material as well as its density and the line profile function 
of the absorption coefficient. Since the absorbers are mov- 
ing in relativistic speeds, the bulk velocity is an important 
factor. Lorentz contraction increases the absorption coef- 
ficient accordingly, and Doppler shift alters the frequency 
of the emission as seen by the absorbers. 

The potential energy liberated by material infalling 
into a black hole is of order a few percent of the material's 
rest mass, and the energy corresponding to 'thermal' kine- 
matic velocity dispersion approaches this rest mass energy. 
Thus, we replace the conventional Maxwellian distribu- 
tion by the Jiittner distribution for relativistic particles 
in deriving the resonant line absorption coefficient. This 
distribution does not give a Gaussian absorption profile 
even in the bulk rest frame of the absorbers. 

The absorption line profiles of the direct images show 
a dip around the fine rest frequency (energy) {E/Eq = !)• 
Absorption can change the line profile significantly. The 



flux around the rest frequencies is, however, augmented 
by the flux from shifted lines from the higher order im- 
ages, especially at high disk inclination. These two effects 
compete, and when the line-of-sight optical depth is high, 
the contribution of the high-order images is masked. 

8. Conclusion 

We present a numerical ray-tracing method for radiative- 
transfer calculations in curved space time and apply the 
method to calculate line emission from accretion disks 
and tori around black holes. Our calculations have shown 
that lines from relativistic accretion tori have profiles very 
different to lines from relativistic thin planar accretion 
disks for the same system parameters, such as the spin 
of the black hole and the viewing inclination. The self- 
obscuration of the inner region of the accretion torus leads 
to weaker red wing in the emission lines when compared 
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with the hnes emitted from a thin planar accretion disk. 
At high inchnation angles the strong blue peak is also ab- 
sent in the line emission from the tori. 

We also investigate the effects of resonant absorp- 
tion/scattering by the line-of-sight material in relativistic 
motion with respect to the emitters in the disk/torus, and 
the observer. Our method does not invoke the Sobolov ap- 
proximation, and the resonant absorption/scattering con- 
dition is derived directly. We have shown that absorp- 
tion effects are important in shaping the profiles of emis- 
sion lines. The interpretation of observations of relativistic 
lines from AGN is non-trivial when absorption is present. 
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Appendix A: Frequency Bins 

Using an astrophysical model together with the metric, we 
can obtain an equation describing the flow of the medium 
through which the ray is propagating; u°^{x°'{X)). The 
properties and distribution of the medium can be used 
to derive the absorption and emission terms xo(2;"(A)) 
and riQ{x°'{X)). These can then be used to calculate what 
is observed from infinity via equation (|17|l . In practice 
this simple method is altered with a ray-tracing numeri- 
cal algorithm. The radiative transfer equation is not time- 
symmetric, so if one follows rays back from the observer 
to the emitter the above formulation cannot be used. The 
emission frequency and intensity is not known until the 
emitter is reached along the path. 

This problem can be solved using one of two ap- 
proaches. The first is to collate the path from emitter to 
observer, and then to integrate along that path forward 
in time using equation (|17|l . The second approach is to 
collate the optical depth for T for a table of frequencies as 
one travels back in time. This optical depth can then be 
used to calculate the observed intensity, /, for each binned 
frequency. The results for each emission element along the 
path are simply summed. 

Since, in general, the resulting frequencies are binned 
anyway to get a spectrum, we choose the second method. 
It requires less storage, and is a simple extension of the 
grey-absorber case where xo and 770 a-i'e no longer functions 
oi I'. In our treatment, we sum over the observed frequen- 
cies rather than those in the rest frame, avoiding the need 
to transform the distribution of the frequencies along the 
path due to the gravitational redshift factor. The act of 
going from the emission frame to the observer's frame also 
affects the intensity /, so we use the invariant I instead, 
only converting to specific intensity just before outputting 
the results. 

Binning the optical depth with frequency along a path 
also allows us to model the radiative transport of the con- 
tinuum as well as the lines, provided that one is in the 
limit where scattering is unimportant. (Stimulated emis- 
sion can be modeled by using a negative absorption coef- 
ficient.) Depending on the relative intensities of the con- 
tinuum and the line, absorption may cause the emission 
feature to be converted into an absorption feature. 

We assume that the continuum is a power law, which 
we parametrize by a slope and intensity at the line rest 
frequency. We have assumed that the continuum is emit- 
ted with an intensity that scales with the line emissivity. 
This means that the equivalent width of the emission line 
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is constant across the disk. We set this to be 0.05 of the 
rest line frequency. (This vahie is roughly what is seen in 
AGN.) To simplify things further, we will fix the power 
law index of the continuum to be 7 = 0.5, where 7 is 
defined by 



I ^CE- 



(A.l) 



in which / is the continuum intensity as a function of en- 
ergy, E^ and C is derived from the given equivalent width. 
This rather hard spectrum was chosen to emphasise the 
line. 

The treatment assumes that the continuum is created 
in a relatively thin planar structure above the accretion 
disk. In effect, we treat the emission from the disk corona 
as part of the injection spectrum, together with the emis- 
sion line, which we propagate through the absorbing ma- 
terial suspended much higher above. 

Appendix B: Ray Tracing Algorithm 

A direct ray-tracing method is used instead of the con- 
ventional transfer-function method, as it is easier to in- 
corporate the numerical radiative-transport calculations. 
The ray-tracing algorithm is as follows: 

1. Integrate the equations governing the geodesies, and 
those describing the optical depth for each frequency, 
from the observer to the emitting surface; 

2. At each crossing of the equatorial plane / torus surface, 
collate the position and the direction of the photon; 

3. Construct the image, and determine the observed fre- 
quency/energy shift; 

4. Integrate the emission over the images of each order to 
produce the line profiles. 

We use a Runge-Kutta integrator to calculate the ray 
trajectories. The relative tolerance is set to 10""'^'^, which 
prevents the ray tracer from missing the torus by taking 
steps which are too large. It also allows a simple Eulerian 
method to be used to integrate the optical depths. The 
absorption coefficient is a slowly varying function on the 
scale of the step size, which allows this optimisation. 

The foot points of the null geodesies (photon trajecto- 
ries) on the disk surface are calculated by a root-finding 
algorithm (see e.g. Press et al. 1992, p. 343). The first four 
intersections of the null geodesies and the disk plane (cor- 
responding to the direct and first three higher-order im- 
ages), and the four- vectors of the photons emitted from 
there are recorded. The incorporation of an emissivity law 
is therefore straightforward, as it is defined in terms of the 
spatial coordinates on the disk plane. Since the trajectory 
of the photon is also saved at each crossing point, it is also 
possible to include the effects of limb darkening, and to 
model a semi-transparent disk. 

Since the disk is imaged upon a sky plane, all grav- 
itational lensing effects on the intensity of the light are 
implicitly included in the calculation of the image itself. 
If a region of the disk is magnified then it will cover more 



area in the image, and thus will appear brighter than a 
non-magnified region. Inclination effects are also included 
implicitly. An inclined disk will cover less pixels than a 
face-on disk, where the number of pixels is roughly pro- 
portional to cosi, where i is the inclination angle. (Light 
bending causes this Euclidean formula to be only an ap- 
proximation.) Gravitational lensing does not alter the ob- 
served surface brightness of a point. 

This implicit inclusion of changing areas of photon 
flux-tubes linking the observer to the emitter vastly sim- 
plifies the calculation of the observed flux. All that is re- 
quired is to integrate over each pixel on the image, taking 
into account the redshift of the emission regions corre- 
sponding to the pixels. If one were integrating over the 
surface of the disk, instead of over the image, then the 
Jacobian of the transformation from the disk to the image 
plane coordinates would be required. This is numerically 
difficult to obtain, and would require a separate transfor- 
mation for each image order. 

The programs (written in Fortran) take a few hundred 
kilobytes of memory under Tru64 Unix on a OSFl V5.1 
1885 Alpha. One version outputs the intensity / frequency 
shift of the emission line for each 'pixel' in the sky plane. 
Another version collates the continuum as well as the line, 
and outputs a table of intensity verses energy integrated 
over the whole disk. Depending on load on the shared 
system, and on the image size requested, a typical run 
takes a few hours per image. 

Appendix C: Surface Finding Algorithm 

We consider the following algorithm to determine the 
torus surface. We integrate equation if^ and tabulate 
the resulting points (r, 0) along the path of the integra- 
tion. Then we interpolate (r, 9) and construct the torus 
surface, where the emission originates. We use spline in- 
terpolation between the surface points. 

When ray tracing the photon paths, we determine the 
intersection of the trajectory and the torus surface. As 
the photon trajectory calculations may take large spatial 
steps, there is a possibity that the torus is not 'detected'. 
To prevent this from happening, we consider the following 
procedure. We take note of the region where the photon 
is located during the trajectory calculation: either inside 
or outside the torus and either above or below the equa- 
toral plane. Whenever the photon leaves one of the four 
regions, and enters another region, we use a boundary- 
searching algorithm to find the exact location where the 
transit occurs. If the trajectory hits the equatorial plane 
outside the torus, the integration will continue. If the tra- 
jectory hits the torus, then integration is terminated. 

By taking smaller steps required by the boundary- 
finding algorithm, we can prevent the integrator from 
missing the torus entirely. This algorithm can be used for 
more complicated surfaces. It works well, because the in- 
tegrator will only miss intersections when the trajectory 
is close to tagential to a surface. This happens close to 
the equatorial plane in the torus models. Adding in a fake 
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boundary there, and thus decreasing the step size, helps 
in preventing missed intersections. 



